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Abstract 

We present the full set of power spectra of cosmic microwave back- 
ground (CMB) temperature and polarization anisotropies due to the 
coupling between quintessence and pseudoscalar of electromagnetism. 
This coupling induces the rotation of the polarization plane of CMB 
photons, thus resulting in non-vanishing B mode and parity-violating 
TB and EB modes. Using the 2003 flight of BOOMERANG (B03) 
data, we give the most stringent constraint on the coupling strength. 
In some cases, the rotation-induced B mode can confuse the hunting 
for the gravitational lensing-induced B mode. 

The existence of a dark component with an effective negative pressure, 
supported by several observations especially the Hubble diagram for the 
type-la supernovae (see, for example, [1, 2]), is still one of the puzzles in 
cosmology. Cosmological constant is the simplest possibility for such dark 
component. However, the observed value of the cosmological constant is 
completely different from theoretical expectation [3]. An alternative candi- 
date, the so-called quintessence described by a dynamical scalar field cf>, is 
naturally considered. The dynamics of (j) in general quintessence models is 
governed by a scalar potential V{4>) which makes the dark energy dominant 
in recent epoch. There are many different kinds of proposed potentials, for 
example, pseudo Nambu-Goldston boson, inverse power law, exponential, 
hyperbolic cosine, and tracking oscillating [4]. To differentiate between the 
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models and finally reconstruct V{<f>) would likely require next-generation 
observations. 

The quintessential potential V{4>) and the field (f> itself are difficult to be 
measured directly. What we can do is to investigate the dark energy density 
0,0 and the time evolution for the equation of state (EOS) = P^/p^, both 
of which are governed by the dynamics of 4>. Several observations have been 
made to give the answers: 157 supernovae in a redshift interval, 0.015 < z < 
1.6, in the "Gold Sample" obtained from a combination of ground-based 
data and the Hubble Space Telescope [2] as well as 115 supernovae with 
0.015 < z < 1 from the Supernova Legacy Survey provide information about 
the dark energy [5] . It is very difficult to determine whether the quintessence 
is more preferred than the cosmological constant by observations of the CMB 
temperature anisotropy spectrum only. Joint analysis of CMB data with 
supernovae or/and large scale structure survey such as SDSS or 2dfGRS 
can offer better constraints on quintessence models [6]. Recently, the study 
of the cross-correlation of maps between CMB and various tracers of matter 
through the integrated Sachs- Wolfe effect has also been carried out by several 
groups [7]. However, the intrinsic properties of dark energy are not well 
constrained so far from the investigations above, thus allowing a very wild 
range of the EOS, which is strongly model dependent. 

An alternative way to study the quintessence is to consider the interac- 
tion of (f) to ordinary matter. Coupling of 4> to dark matter is considered as a 
possible solution for the late time coincidence problem [8]. It leaves distinct 
imprints on the CMB temperature anisotropy and the matter power spec- 
trum due to an excess of cold dark matter at early epoch when compared 
to the standard cosmology model [9] (We define the "standard cosmology" 
as the model without coupling between quintessence and ordinary matter). 
Of particular interest in this Letter, we study the coupling of (f> to the pseu- 
doscalar of electromagnetism. Pseudoscalar couplings usually arise from the 
spontaneous breaking of a compact symmetry group, say, U(l) (see Erie- 
man et al. 1995 at reference [4]). Carroll has argued that the coupling 
P F p4>/ ' MF^ P F^ V leads to the rotation of the polarization vector of propa- 
gating photons if 4> is varying with time. Here /3 F p is the coupling strength, 
M is the reduced Planck mass, and F^ v is the dual of the electromagnetic 
tensor. This effect is called as the "cosmological birefringence" [10]. 

Measurements of the polarization of distant astronomical objects would 
provide information about the cosmological birefringence. Carroll used the 
rotation of polarization direction for distant radio sources to constrain the 
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coupling strength [10]. Another proposed method is the CMB polarization 
from the last scattering surface [11, 12]. In this Letter, we study the effects 
of the cosmological birefringence on the power spectra of CMB temperature 
and polarization. Only assuming the shape of the potential V(4>), we present 
the first CMB power spectra in the presence of the cosmological birefringence 
by the use of the full Boltzmann code. 

Thomson scatterings of anisotropic radiation by free electrons give rise to 
the linear polarization, which is usually described by the Stokes parameters 
Q and U [13]. In standard cosmology, the time evolution of the polarization 
perturbation is governed by the Boltzmann equation [14]. If there is a physi- 
cal mechanism which rotates the polarization plane, the evolution equations 
for the Fourier modes of the Stokes parameters are modified to 



AQ±iu(Kv) + iknA Q ± iU Qs.,r)) = n e a T a{rf) 

na)4 m) (M) 




-AQ±j[/(k, rj) 
Ti2ujA QTiU (k,ri), (1) 



where the derivatives are taken with respect to the conformal time rj, /j, = nk 
is the cosine of the angle between the CMB photon direction and the Fourier 
wave vector, n e is the number density of free electrons, ot is the Thomson 
cross section, and a is the scale factor. s YJ m is spherical harmonics with 
spin-weight s with m = 0, ±1, ±2 corresponding to scalar, vector, and tensor 
perturbations, respectively, if the axis of s YJ m is aligned with the wave vector 
k. 4 m) is the source term of generating polarization, which is the compo- 
sition of the quadrupole components of the temperature and polarization 
perturbations S l ™\\s,r)) = A^ } (k,r/) + 12>/6Aj$ (k, rj) + 12>/6A^(k, rj). 
We have followed the notation in Ref. [15] and expanded the perturbations 
in the spin-0 and spin-2 spherical harmonics [16] according to the scalar 
and tensor properties of the temperature anisotropy and polarization re- 
spectively. Thus, Aj^ and A^ are the respective expansion coefficients 
for Ay on basis Yj jTn and A.Q±m on basis ±riXim- The last term in Eq. (1) 
appears due to the rotation of the polarization plane. The dispersion rela- 
tion for electromagnetic radiation coupling to the time varying quintessence 
field (p is given by E 2 = k 2 ± k(3 F p(j)/(aM), where ± refer to the right and 
left handed circular polarization, respectively. Therefore, the net angular 
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velocity of the polarization plane is [10] 

We are used to decomposing the polarization on the sky into a divergence 
free component, the so-called E mode, and a curl component, the so-called 
B mode because the values for Q and U depend on the choice of a coordinate 
system. Whether B mode is generated depends on the existence of U for 
the local mode whose wave vector k parallels to the z of the coordinates 
whereas Q is defined as the difference in intensity polarized in the 6 and 
ip directions [17]. Mathematically, for m = 0, only Q is generated in the 
local mode. The axisymmetry of the radiation field about the mode axis 
guarantees that no B mode can be generated by scalar mode perturbations. 

In the presence of cosmological birefringence, we can find two important 
features in Eq. (1). Firstly, the rotation of the polarization plane generates 
U contributions to the local mode polarization. This converts the power 
from the E mode to B mode. The conversion depends on how much rota- 
tion there is from the epoch when the polarization is generated to today. 
Secondly, TB and EB cross correlations are expected to vanish due to the 
parity (T and E have parity (— 1) while B has (— 1)^ +1 ). The cosmological 
birefringence violates the parity and thus generates the TB and EB power 
spectra whose magnitudes depend on the integrated rotation of polarization 
too. Substituting the angular velocity of the polarization plane in Eq. (2) 
into Eq. (1), we calculate the power spectra T, E, B, TE, TB, and EB 
modes. 

These six power spectra form a complete two-point statistics of CMB 
temperature and polarization anisotropics. To simplify the calculation, we 
only include the scalar perturbations. That is by setting m = in Eq. (1). 
Without showing the details, we just write down the power spectra which 
are obtained from the solutions for the line-of-sight integration as 
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Cj E = (4nfJ^^^Jk'dkA T (k,r,o)A E (k, Vo ), 



Cj B = (4nfJ^^±J-Jk 2 dkA T (k, V o)A B (k, Vo ), (3) 
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where 



A E (k,r] ) +iA B (k,rjo) = / drjg(r))S P (k,r]) e 
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where the visibility function g(r]) describes the probability that a photon 
scattered at epoch r] reaches the observer at the present time, ijq. Similar 
to Sp = Sp^ , St is the source term generating the temperature anisotropy. 
ji is the spherical Bessel function and r = r/Q — -q. The rotation angle 
a(r]) = dr]'uj(rj'). We do not present the formula for the temperature 
anisotropy because it is unchanged under the rotation of the polarization 
plane. 

We have modified the public code CMBFast [18] for our purpose. Here, 
we consider the potential V((f>) = Vq exp(A0 2 /2M 2 ) for our quintessence 
model (the hyperbolic cosine potential is also considered, see below), where 
A is a parameter determining how shallow the potential is. Hereafter we 
fix A = 5 and we will obtain similar results by choosing other values of A. 
We plug the table of the EOS for this quintessence model into the modified 
CMBFast code and input the cosmological parameters from the best fit 
values of the WMAP three-year results [19]. The power spectra are then 
normalized to the first peak of the temperature anisotropy measured by the 
three-year WMAP observation [20]. The left panel of Fig. 1 shows the E 
and B mode power spectra with the coupling strength (3 F p ranging from 
10~ 5 to 1(T 3 . The EB mode power spectrum is shown in the right panel. 
On small scales, increasing coupling strength results in a suppression of 
the E mode in the standard model and non- vanishing B and EB modes. 
Furthermore, the shapes of the B and EB mode power spectra basically 
follow the standard E mode except the reionization bump on large scales. To 
explain this, we may make a rough estimation in Eq. (3): Cbi ~ Ce£ sin 2 2ot£ 
and Ceb£ ~ 0-5Ce£ sin4c^ with ot£ being the total rotated angle for certain 
angular scales 9 ~ tt/£ from last scattering epoch to today. This ct£, in 
general, is not constant for all the scales. The E mode power on small scales 
mainly comes from the recombination epoch at z ~ 1100. On the other hand, 
the boosting power on large scales comes from reionization epoch when the 
CMB photons are rescattered by free electrons at z ~ 10 [21]. From Eq. (2) 
and the evolution of (f>, we find that the integrated rotation angle from the 
reionization epoch is much smaller than that from the recombination epoch. 
Therefore, there is much less power converted from E mode to B mode on 
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large scales than small scales. 




Figure 1: E, B (left panel; the lower three thick curves are B modes) and 
EB (right panel) mode power spectra from the cosmological birefringence 
with different coupling strength. 

If the coupling strength is large enough, the B mode induced by the 
cosmological birefringence will mix up with the gravitational lensing-induced 
B mode. Gravitational lensing by large scale structures modifies slightly the 
primary E mode power spectrum. Most noticeably it generates, through 
mode coupling, B mode polarization out of pure E mode signal [22]. The 
lensing-induced B power spectrum, which peaks around I ~ 1000, has the 
roughly similar shape with that from the birefringence. We also show the 
power spectrum of the lensing-induced B mode in Fig. 1 by a thin solid curve 
for comparison. The birefringence-induced B mode is indeed compatible 
with the lensing-induced B mode for (5 F p ~ 10 -4 . 
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Figure 2: TE (left panel) and TB (right panel) mode power spectra from 
the cosmological birefringence with different coupling strength. 

Fig. 2 shows the TE and TB power spectra for different coupling strength. 
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Having the complete set of power spectra, we can constrain the coupling 
strength from observed data. In order to focus on the cosmological birefrin- 
gence, we do not make the global fit to all the cosmological parameters. It 
is debatable that f3 p p is degenerate with other cosmological parameters due 
to the decrement of the TE and E power spectra. As we will see later, the 
upper limit on the coupling strength prevents it from making a significant 
change on the TE and E modes. 




Figure 3: Likelihood function of the coupling strength. 

We use the North American data from the 2003 Antarctic flight of 
BOOMERANG [23] to calculate the X 2 = - T i) C iJ{ D ) ~ T j) fit- 

ting for (5 F p while fixing other cosmological parameters, where D\ is the 
ith band-power, C is the covariance matrix, and Zf = Ylt^tWit/^ is ob- 
tained by the theoretical prediction multiplied by the band-power window 
function Wajl. Those band-power data, including covariance matrices and 
window functions, are available; see [24]. In practice, we only fit (3 FF to the 
EB and TB power spectra which come from the parity violation. We do 
not use the B mode power spectrum because it is contaminated by lensing- 
induced signal. If the lensing-induced B mode can be successfully cleaned 
by appropriate techniques such as those proposed by Seljak and Hirata [25], 
we expect that including the B mode will give a stronger constraint. We 
convert % 2 to the likelihood by C = e~ x I 2 and normalize the maximum 
likelihood value to unity. The result is shown in Fig. 3. The upper limit on 
the coupling strength is found to be \P FF A(f>\/M < 8.32 x 10 -4 at 95% con- 
fidence level, where A<p is the total change of 4> till today. This small value 
of the coupling strength gives an insignificant change on TE and E modes 
and thus will not affect the determination of the cosmological parameters. 
We also use the hyperbolic cosine potential V(4>) = Vq cosh(A</>/M) for the 
testing. Even though the quintessence evolution in this potential is different 
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from that in the exponential case, the upper limit value of \f3 FF A<fi\/ M does 
not change much. This value is much smaller than the result in Ref. [10], 
3 x 10~ 2 , where Acj) is only from z = 0.425 to today. It is remarkable that 
the rotation-induced B mode with the upper limit value of the coupling 
strength exceeds the lensing-induced B mode. Therefore, careful measure- 
ments of TB and EB are necessary for separating the two effects. 

Several authors have studied the effect of parity violation on the CMB 
power spectra by assuming a constant rotation angle a [11, 12]. They ob- 
tained a new set of rotated power spectra from combining the power spec- 
tra in the standard model with the sine or cosine function of a. Further- 
more, taking the rotation angle as a free parameter, Feng et al. [12] con- 
strained the rotation angle using the data made by the first year WMAP 
and BOOMERANG observations. However, the time-varying scalar field <j> 
in their work is constrained such that the integrated rotation angle should 
be very small from the recombination to the reionization epoch. Therefore, 
it is less supported by general quintessence models. 
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